Flavor coupled with chiral oscillations in the presence of an 

external magnetic field 

A. E. Bernardini 

Institute) de Fisica Gleb Wataghm, UNICAMP, 
PO Box 6165, 13083-970, Campinas, SP, Brasil* 
(Dated: February 1, 2008) 

Abstract 

By reporting to the Dirac wave-packet prescription where it is formally assumed the fermionic 
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introduce some small modifications to the standard neutrino flavor conversion formula. Assuming 
the correspondent spinorial solutions of a "modified" Dirac equation, we are specifically interested 
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packet treatment becomes clear when we assume B orthogonal/parallel to the direction of the 
propagating particle. 
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I. INTRODUCTION 



Obtaining exact solutions of a generic class of Dirac wave equations [1-5] becomes im- 
portant since, for many times, the conceptual understanding of physics can only be brought 
about by such solutions. These solutions also correspond to valuable means for checking and 
improving models and numerical methods for solving complicated physical problems. In the 
context in which we intend to explore the Dirac formalism, we can report about the Dirac 
wave-packet treatment which can be useful in keeping clear many of the conceptual aspects 
of quantum oscillation phenomena that naturally arise in a relativistic spin one-half particle 
theory. These studies have been paralleled by much progress on the theoretical front of the 
quantum mechanics of neutrino oscillations, both in phenomenological pursuit of a more 
refined flavor conversion formula [6-8] which, sometimes, deserves a special attention, and 
in efforts to give the theory a formal structure within quantum field formalism [9-11]. Un- 
der the point of view of a first quantized theory, the flavor oscillation phenomena discussed 
in terms of the intermediate wave-packet approach [12, 13] eliminates the most controver- 
sial points rising up with the standard plane-wave formalism [14, 15]. However, a common 
argument against the intermediate wave-packet formalism is that oscillating neutrinos are 
neither prepared nor observed [16]. This point was clarified by Giunti [10] who suggested a 
solution by proposing an improved version of the intermediate wave-packet model where the 
wave packet of the oscillating particle is explicitly computed with field-theoretical methods 
in terms of external wave-packets. Since we intend to concentrate the discussion on the Dirac 
equation properties for rigorously deriving a flavor/chiral conversion formula for fermionic 
particles non-minimally coupling with an external magnetic field, in this preliminary inves- 
tigation, we avoid the field theoretical methods in detriment to a clearer treatment with 
intermediate wave-packets which commonly simplifies the understanding of physical aspects 
going with the oscillation phenomena. 

From the experimental point of view, compelling evidences have continuously ratified 
that neutrinos undergo flavor oscillations in vacuum or in matter. For instance, we refer to 
the outstanding results of the Super-Kamiokande atmospheric neutrino experiment [17], in 
which a significant up-down asymmetry of the high-energy muon events was observed. We 
also have the results of the SNO solar neutrino experiment [18, 19] in which a direct evidence 
for the transition of the solar electron neutrinos into other flavors was obtained. Finally, the 
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KamLAND experiment [20] has confirmed that the disappearance of solar electron neutrinos 
is mainly due to neutrino oscillations and not to other types of neutrino conversions [21, 
22]. The experimental data could be completely interpreted and understood in terms of 
three flavor quantum numbers, excepting by the LSND anomaly [19, 23, 24] which permit 
us to speculate the existence of (at least) a fourth neutrino (flavor?) which has to be 
inert. In fact, the hypothesis of mixing between known neutrino species (electron, muon 
and tau) and higher mass neutrinos including sterile neutrinos was naively studied in the 
literature [25, 26]. The neutrino spin-flipping attributed to some dynamic external [27] 
interacting process, which comes from the non-minimal coupling of a magnetic moment with 
an external electromagnetic field [28], was formerly supposed to be a relevant effect in the 
context of the solar-neutrino puzzle. As a consequence of a non- vanishing magnetic moment 
interacting with an external electromagnetic field, left-handed neutrinos could change their 
helicity (to right-handed) [29] . The effects on flavor oscillations due to external magnetic 
interactions in a kind of chirality-preserving phenomenon were also studied [30] but they 
lack of a full detailed theoretical analysis. In the standard model flavor-changing interactions, 
neutrinos with positive chirality are decoupled from the neutrino absorbing charged weak 
currents [31]. Consequently, such positive chirality neutrinos become sterile with respect 
to weak interactions. Independently of any external electromagnetic field, since neutrinos 
are detected essentially via V-A charged weak currents, the chiral oscillation mechanism by 
itself may even explain the "missing" LSND data. Despite the experimental circumstances 
not being completely favorable to such an interpretation, which should be an additional 
motivation to our theoretical calculations, the quantum transitions that produces a final 
flavor eigenstate corresponding to an active-sterile quantum mixing is perfectly acceptable 
from the theoretical point of view. 

To summarize, the main point to be considered in this manuscript concerns with accu- 
rately obtaining the complete flavor/chiral conversion formula for fermionic particles, inde- 
pendently of its range of applicability which certainly shall be determined in terms of more 
accurate experimental resolutions. After such an introductory perspective, the first step of 
our study, which is presented in section II, concerns with the immediate description of chiral 
oscillations in terms of the Hamiltonian dynamics ruled by the Lagrangian for the inter- 
action between a fermionic field ip{x) and an electromagnetic field F^ v (x). We investigate 
how the interactions with an external magnetic field can change the characteristics of chiral 
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oscillations which were previously obtained for free propagating particles in vacuum [8]. By 
including the effect of chiral oscillation, our final aim consists in verifying how the interac- 
tion with an external magnetic field can modify the neutrino flavor oscillation formula which 
was previously established by means of what we call scalar prescription [32]. In section III, 
for treating the time-evolution of the spinorial flavor-eigenstates, we shall take into account 
the chiral nature of charged weak currents and the time-evolution of the chiral operator 
with Dirac wave-packets. To do it, we shall use the "modified" Dirac equation (with the 
above non-minimal coupling term) as the evolution equation for the mass-eigenstates. We 
just remark that we do not mind in developing the calculations concerning with analytic 
characteristics of the localization of each mass(flavor)-eigenstate since all the modifications 
introduced by this prescription can be read independently of the wave-packet shapes. How- 
ever, we also dedicate an appendix to a small revision of the fundamental points of the 
wave-packet formalism which, eventually, can be important for proceeding with subsequent 
analytic calculations. Finally, we draw our conclusions in Section IV. 

II. CHIRAL OSCILLATIONS IN THE PRESENCE OF AN EXTERNAL MAG- 
NETIC FIELD 

In order to introduce the coupling with external magnetic fields, we observe that even 
presenting an electric charge neutrality, neutrinos can interact with a photon through loop 
(radiative) diagrams. The Lagrangian for the interaction between a fermionic field ip(x) and 
an electromagnetic field written in terms of the field-strength tensor F^ u (x) = d tx A v {x) — 
d v A tx (x) is given by 

C = ^#r) [<v (fxF^(x) - d^{x)) + h.c] ij>(x) (1) 

where x = (t,x), a^ v = |[7 M ,7i/], the dual field-strength tensor T^ v {x) is given by 
T^ v (x) = | e ^ A<5 F xs (x)[43] and the coefficients /z and d represent, respectively, the mag- 
netic and the electric dipole moment which establish the neutrino electromagnetic coupling. 
One can notice that we have not discriminated the flavor/mass mixing elements in the above 
interacting Lagrangian since we are indeed interested in the physical observable dynamics 



4 



ruled by the Hamiltonian 

H = a-p + pm-p (// F^{x) - + h.c. 

= a-p + ^m-pS- B(x) - d £ • £(z)] , (2) 

where, in terms of the Dirac matrices, a = X)k=i a kk = YH=i lolkk, (3 = 7 , and and 
E{x) are respectively the magnetic and electric fields. In fact, the Eq. (2) could be extended 
to an equivalent matrix representation with flavor and mass mixing elements where the 
diagonal (off-diagonal) elements described by fj, itj (m itj ) and d itj (m itj ), where i, j are mass 
indices, would be called diagonal (transition) moments. In this context, for both Dirac and 
Majorana neutrinos, we could have transition amplitudes with no n- vanishing magnetic and 
electric dipole moments [33-35]. Otherwise, the CP invariance holds the diagonal electric 
dipole moments null [36]. Specifically for Majorana neutrinos, it can be demonstrated that 
the diagonal magnetic and electric dipole moments vanish if CPT invariance is assumed 
[33]. 

Turning back to the simplifying example of diagonal moments and assuming CP and CPT 
invariance, we can restrict our analysis to the coupling with only an external magnetic field 
B(x) by setting d = 0. From this point, the expression for jj, can be found from Feynman 
diagrams for magnetic moment corrections [36] and turns out to be proportional to the 
neutrino mass (matrix), 

3eG 3m e G n _ ,„_ in m v . . 

a = — — — m = — = — u s m v = 2.7 x 10 a B (3) 

8V2n' 2 4V2n 2 m N 

where G is the Fermi constant and m N is the nucleon mass [44]. In particular, for m u m 1 eV, 
the magnetic moment introduced by the above formula is exceedingly small to be detected 
or to affect astrophysical or physical processes. 

Since we are interested in constructing the dynamics ruled by the Hamiltonian of Eq. (2), 
we firstly observe that the free propagating momentum is not a conserved quantity, 

j f (p) = i([H,p}) = (£•£?(,))> (4) 

In the same way, the particle velocity given by 

j t (x) = i([H,x}) = (a) (5) 

comes out as a non-null value. After solving the x(t)(a(t)) differential equation, it is possible 
to observe that, in addition to a uniform motion, the fermionic particle executes very rapid 



oscillations known as zitterbewegung [38]. By following an analogous procedure for the Dirac 
chiral operator 75 , newly recurring to the equation of the motion, it is possible to have the 
chirality and the helicity dynamics respectively given by 

^( 7 5 > = 2*( 7o75 [m- /U S-B(x)]> (6) 

and 

j f (h) = ^(7oP'V)(S-B(,)) + 2(SxB(,)).p]> (7) 

where we have alternatively defined the particle helicity as the projection of the spin angular 
momentum onto the vector momentum, h = |S • p (with p in place of p). From Eqs. (6-7) 
we can state that if a neutrino has an intrinsic magnetic moment and passes through a region 
filled by an magnetic field, the neutrino helicity can flip in a completely different way from 
how chiral oscillations evolve in time. In the non-interacting case, it is possible to verify 
that the time-dependent averaged value of the Dirac chiral operator 7 5 has an oscillating 
behavior [31] very similar to the rapid oscillations of the position. The Eqs. (6-7) can be 
reduced to the non- interacting case by setting B(x) — so that 

j t (h) = i([H,h}) = -((axp).p) = (8) 

and 

j t (f) = z([H,r']) = 2im<7V>. (9) 

from which we confirm that the chiral operator 7 s is not a constant of the motion [31]. 
The effective value of Eq. (9) appears only when both positive and negative frequencies are 
taken into account to compose a Dirac wave-packet, i. e. the non-null expectation value of 
(7o7s) is revealed by the interference between Dirac equation solutions with opposite sign 
frequencies. The effective contribution due to this interference effect lead us to report to the 
Dirac wave-packet formalism in order to quantify neutrino chiral oscillations in the presence 
of an external magnetic field. 

Assuming the simplifying hypothesis of a uniform magnetic field B, the physical impli- 
cations of the non-minimal coupling with an external magnetic field can then be studied by 
means of the eigenvalue problem expressed by the Hamiltonian equation 

H( P ) <p n = E n ( P ) (p n — {oc-p + P[m- /j,~E- B]}ip n , (10) 
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for which the explicit 4x4 matrix representation is given by 



H(p)<p n 



m - fiB t -fi(B x - iB y ) p z p x - ip y 

-H{B X + iB y ) m + fiB z p x + ip y -p z 

p z Px-iPy -(m-fj,B z ) /i(B x -iB y ) 

p x + ip y -p z n(B x + iB y ) -{m + nB x ) 



(11) 



The most general eigenvalue (E n (p)) solution of the above problem is given by 

E n (p) = ±E S ( P ) = ±^m 2 + p 2 + a 2 + (-iy2^m*a> + (p x a7, s = 1, 2, (12) 

where we have denoted E n=1 ^ 3A = ±E 3=lt2 and we have set a = fiB. The complete set 
of orthonormal eigenstates <p„ thus can be written in terms of the eigenfunctions U{p s ) 
with positive energy eigenvalues (+E s (p)) and the eigenfunctions V(p s ) with negative energy 
eigenvalues (—E s (p)) : 



U( Ps ) = N( Ps ) 
V(p.) = N( Ps ) 




at I A 7 at 




'A- a- 

s s 

A+a+' 



(13) 



where p s is the relativistic quadrimomentum, p s = (E s (p),p), N(p s ) is the normalization 
constant and 

A± = A» ± 2m|a| - a 2 , = 2E s (p)\a\ ± (A 2 ( P ) + a 2 ) 



with 



A 2 (p) = 2 (a 2 + (-l)V«i 2 a 2 + (pxa) 2 ) . 



(14) 



We can observe that the above spinorial solutions are free of any additional constraint, 
namely, at a given time t, they are independent functions of p and they do not represent 
chirality/helicity eigenstates. 

In order to describe the above Hamiltonian dynamics for a generic observable 0(t) we 
can firstly seek a generic plane-wave decomposition as 



exp [—i(E 3 (p) t — p ■ x)] U (ps), for positive frequencies and 
exp [ i(E s (p) t — p • x)} V(ps), for negative frequencies, 



(15) 
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so that the time-evolution of a plane-wave-packet ip(t,x) can be written as 

<$p 



if;(t, x) 



(2tt)' 



J2{b(Ps)U(Ps) exp \-iE s { P )t] 

s=l,2 

+ c?*(p s )V(p s ) exp [+i £ s (p) t}} exp [ip • cc], 



(16) 



with p s = (E s , —p). The Eq. (16) requires some extensive mathematical manipulations for 
explicitly constructing the dynamics of an operator 0{t) like 



0(t) 



X 



X 



(17) 



If, however, the quoted observables like the chirality 7 5 , the helicity h or even the spin 
projection onto B commuted with the Hamiltonian H, we could reconfigure the above 
solutions to simpler ones. To illustrate this point, let us limit our analysis to very restrictive 
spatial configurations of B so that, as a first attempt, we can calculate the observable 
expectation values which appear in Eq. (6). Let us then assume that the magnetic field 
B is either orthogonal or parallel to the momentum p. For both of these cases the spinor 
eigenstates can then be decomposed into orthonormal bi-spinors as 



Uip a ) = N + ( Pa 



and 



V( Pa ) = N~(p s 



_X + (Ps) 

X~(Ps) 



(18) 



(19) 



Eventually, in order to simplify some subsequent calculations involving chiral oscillations, 
we could set <^* 2 and Xi, 2 as eigenstates of the spin projection operator <x • B, i. e. beside 
of being energy eigenstates, the general solutions U(p s ) and V(p B ) would become eigenstates 
of the operator S • B and, equivalently, of X • a. 

Now the Eq.(lO) can be decomposed into a pair of coupled equations like 



± 

s 1 



(±E S - m + a ■ a) ip* = ±cr p\ 
(±E S + m - cr ■ a) xf = ±tr-p<pf, 



(20) 



where we have suppressed the p s dependence. By introducing the commuting relation [cr ■ 
p, a ■ B] =0 which is derived when p x B = 0, the eigenspinor representation can be 
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reduced to 



and 



U{ Pe ) 



V{ Pa ) 



E s + m s 



2E. 



E s + m s 



2E« 



E s +m s r a 



(TO 

— Y 

Es+m s A s 



(21) 



(22) 



with m s = m — (-l) s |a| and the energy eigenvalues 



±E S = ±^p T T 



(23) 



In this case, the closure relations can be constructed in terms of 



3=1,2 s=l,2 s 



1 - (-i) s E • a 



1 - (-l) a E • a 



2£„ 



(24) 



s=l,2 s=l,2 

Analogously, by introducing the anti-commuting relation {<r ■ p, cr ■ B} when p ■ B = 0, the 
eigenspinor representation can be reduced to 



U{ Ps ) 



and 



V(p.) = 



e + m 
2e 



£ + m 



2£ 



cr p , 

— <P 

so+m • a 



crp 

— — y 

so+m 



(25) 



(26) 



with e = ^/p 2 + m 2 and the energy eigenvalues 



±E. = ±[e + (-l)'\a\]. 



In this case, the closure relations can be constructed in terms of 



(27) 



53W(p.)®ZAp.)7o 

8 = 1,2 

- J] V(p.) <g> V f (p s )7o 

8 = 1,2 

where p = (e ,p)- 



7^Pu + m 

&= 

~7mPq + m 
2£ 



E 

1,2 

E 

8 = 1,2 



1 - (-l) s 7o E • a 



1 - (-l) s 7 () E • a 



(28) 
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Since we can set y+ 2 = Xi, 2 as ^ ne components of an orthonormal basis, we can immedi- 
ately deduce the orthogonality relations 

U\p s )U(p r ) = VHps)V(p r ) = 5 sr , 

U\vs) To V( Pr ) = VHps) iMM = (29) 

which are valid for both of the above cases. 

Finally, the calculation of the expectation value of 7 s (t) is substantially simplified when 
we substitute the above closure relations into the wave-packet expression of Eq. (16). To 
clarify this point, we suppose the initial condition over ip(t,x) can be set in terms of the 
Fourier transform of the weight function written as 

if( P - Pi )w= i b (P*)U(Ps) + d*(Ps)V(p s )} (30) 

8 = 1,2 

so that 

^(0,aj) = J j^y<P{P-Pi)ex-V[ip- x]w (31) 

where w is some fixed normalized spinor. By using the ort honor mality properties of Eq. (29), 
we find [39] 

6(p.) = <p(p-Pi)U*(p,)w, 
d*(p s ) = <p(p- Pi )V*(p.)w. (32) 

For any initial state ip(0,x) given by Eq. (31), the negative frequency solution coefficient 
d*(p s ) necessarily provides a non-null contribution to the time-evolving wave-packet. This 
obliges us to take the complete set of Dirac equation solutions to construct a complete 
and correct wave-packet solution. Only if we consider the initial spinor w being a positive 
energy (E s (p)) and momentum p eigenstate, the contribution due to the negative frequency 
solutions d*(p s ) will become null and we will have a simple expression for the time-evolution 
of any physical observable. By substituting the closure relations of Eqs. (24) and (28) into 
the time-evolution expression for the above wave-packet, the Eq.(16) can thus be rewritten 
as 



vfj(t,x) 



f^yMP- Pi) ex P l*P ■ X \ E { [ C ° S [Es t] ~ 1 ^ 8111 [Es t] \ ( 1 " ( " 1 2 3S d ) } W ^ 

^ ' s=l,2 
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for the first case where E s is given by Eq. (23) and H s = a. ■ p + 7o ra s , or as 

lK*,*) = J^<p(p-P t )e^[ip-x]^{[coB[E.t]-i^ sm[E 3 t]\ (Mzil^o^a) } w(3 4) 

for the second case where e is given by Eq. (27) and H — ct • p + j m. 

Now we can proceed with the calculation of the time-evolution of 75 (t) described by 
Eq. (6) assuming the simplifying hypothesis of a uniform magnetic field B. Once we have 
assumed the neutrino electroweak interactions at the source and detector are (left) chiral 
(-i/;7 M (i — 7 5 )^W M ) only the component with negative chirality contributes to the propaga- 
tion. It was already demonstrated that, in vacuum, chiral oscillations can introduce very 
small modifications to the neutrino conversion formula [8, 31]. The probability of a neutrino 
produced as a negative chiral eigenstate be detected after a time t can be summarized by 

/l X 
d 3 x ^(t,x) -ZJlj,(t,x) = - (1 - ( 75 )(t)) . (35) 

From this integral, it is readily seen that an initial -1 chiral mass-eigenstate will evolve with 
time changing its chirality. By assuming the fermionic particle is created at time t = as a 
-1 chiral eigenstate ( 7s w = -w), in case of [a ■ p, cr ■ B] = (B parallel to p), we can write 

w] J2 { h cos2 [Es t] + tlJ ^r sin [2 Es t] + ik w k sin2 [E ° t] ] ( 1 " ( "T Ed ) } w 

s=l,2 

= (-i)/ -^y(P - p.) E { H *i + ^ sini & *]] (^f 5 ^) »} 

= ^ /(^^-^ElK + Sf^I 2 ^']] »* (i=t¥^) »}, (36) 

^ ' s=l,2 

where we have used the wave-packet expression of Eq. (33) and, in the second passage, we 
have observed that 



w 



7 5 iu = -1, w^[H s , %}w = and H s 7s H s = p 2 - m 2 s . (37) 



The above expression can be reduced to a simpler one in the non-interacting case [31]. 
Due to a residual interaction with the external magnetic field B we can also observe chiral 
oscillations in the ultra- relativistic limit. However, from the phenomeno logical point of view, 

2 

the coefficient of the oscillating term goes with which makes chiral oscillations become 
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not relevant for ultra-relativistic neutrinos [8, 40]. As a toy model illustration, by assuming 
a highly peaked momentum distribution centered around a non-relativistic momentum p t <C 
m s , and with the wave-packet effects practically ignored, the chiral conversion formula can 
be simplified and written as 

P(^a,L —> ^a,L~, t) ~ - (f + cos [2mt] cos [2 \a\ t] + sin [2 rat] sin [2 \a\ t]w^S ■ aw) (38) 

where, in this case, all the oscillating terms come from the interference between positive 
and negative frequency solutions which compose the wave-packets. Turning back to the 
case where {cr ■ p, cr ■ B} — (B orthogonal to p), we observe a phenomenologically more 
interesting result. By following a similar procedure with the mathematical manipulations, 
we can write 



where we have used the wave-packet expression of Eq. (34) and, in addition to w^ s w — -1, 
we have also observed that {H , 7 5 } = 27 5 S • p and, subsequently, u^E • P7 S ■ aw = 0. 

2 

Now, in addition to the non-interacting oscillating term ^ cos [2 e 1] , which comes from 
the interference between positive and negative frequency solutions of the Dirac equation, 
we have an extra term which comes from the interference between equal sign frequencies 
and, for very large time scales, can substantially change the oscillating results. In this case, 
it is interesting to observe that the ultra-relativistic limit of Eq.(39) leads to the following 
expression for the chiral conversion formula, 



which, differently from chiral oscillations in vacuum, can be phenomenologically relevant. 
Obviously, the Eq (40) is reproducing the consolidated results already attributed to neutrino 
spin-flipping [36] where, only when we take ultra-relativistic limit (ra = 0), the chirality 
quantum number can be approximated by the helicity quantum number. However, now it 
was accurately derived from the complete formalism with Dirac spinors. 



(7 8 >(t) = / 




(39) 



w 




(40) 
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III. FLAVOR COUPLED TO CHIRAL OSCILLATIONS IN THE PRESENCE OF 
AN EXTERNAL MAGNETIC FIELD 



After obtaining the time-evolution of the spinorial mass-eigenstate wave-packets in the 
presence of an external magnetic field, we have to immediately observe that we cannot arbi- 
trarily overlook chiral oscillations while treating neutrino flavor oscillations. Consequently, 
in order to quantify the interference of chiral oscillations over neutrino flavor conversion 
processes, the chiral nature of charged weak currents and the time-evolution of the chiral 
operator must be aggregated to the flavor oscillation formula. 

Considering that the main aspects of oscillation phenomena can be understood by study- 
ing a two flavor problem, we can introduce this simplifying condition for computing the 
oscillation probabilities. Thus the time-evolution of flavor wave-packets can be described by 
the state vector 

*&(t,x) = ^(t, x) cos 9 v x + i[) 2 (t, x) sin 6 v 2 

= [ip 1 (t, x) cos 2 9 + tp 2 {t, x) sin 2 9} u a + (t, x) — (t, x)] cos 6 sin 6 vp 

= i/j a (t,x;6)v a + ipp(t,x;0)v f 3, (41) 

where u a and Up are flavor-eigenstates, v x and v 2 are mass-eigenstates and 9 is the mixing 
angle. The probability of finding a flavor state Vp at the instant t is equal to the integrated 
squared modulus of the vp coefficient, 

P(v a ^vp;t) = jd 3 x |^(t )a; ^)| 2 = ^{l-DFo(t)}, (42) 

where Dfo(£) represents the mass-eigenstate interference term explicitly computed in terms 
of 

DFO(t) = ^Jd 3 x [iPl(t,x)ip a {t,x) + ^(t.x^^x)] . (43) 

Once we have assumed we have to take into account the chiral conversion character for 
obtaining a complete description of the flavor conversion mechanism, the most complete 
oscillation probability formula must be written as 

P{y^ - vp, L ;t) = ^1 { Dco(t) - DFCo(t) } (44) 

where the results for Dco(t) corresponding exclusively to chiral oscillations for each mass- 
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eigenstate component can be immediately reproduced from Eq. (35) as 



Dco(t) = \J d " x 



= -(l- Ml^±M^ j (45) 

where the averaged values of j 5 can be, for instance, explicitly calculated in terms of the 
results of the Eqs. (36) and (39) and, at the same time, the mixed flavor and chiral oscillation 
term can be given by 



DFCO(t) = i J d 3 x 



ipl(t,x)—-^-ip 2 (t,x) + ^\{t,x) 1 7 V i(*»») 



(46) 



2 Tzy ' ' r2V ' ' 2 
which deserves a more careful calculation. We shall see then how we explicitly construct 
the complete oscillation formula containing both "flavor- appearance" (neutrinos of a flavor 
not present in the original source) and "chiral-disappearance" (neutrinos of wrong chirality) 
probabilities for both of the particular cases discriminated in section II. The following results 
are obtained after some simple but extensive mathematical manipulations where, again, we 
have imposed an initial constraint which establishes that the normalizable mass-eigenstate 
wave functions ^ lt2 (t,x) are created at time t = as a negative chiral eigenstate (wl t2 j 5 w lt 2 = 
-1). All the subsequent calculations do not depend on the gamma matrix representation. In 
correspondence with the first case of the previous section, where the propagating momentum 
p is parallel to the magnetic field B, we can write 

where we have used the correspondence tp 2 (p) = <p(p - Pi)y(p _ P2) ( see the discussion in the 
Appendix) and = \J p 2 + with i — 1,2 corresponding to the mass indices. In fact, 
in the ultra-relativistic limit, and for not very strong magnetic fields (|a| <C \p\), the contri- 
bution due to the very rapid oscillations which come from the interference between negative 
and positive frequency components, analogously to the case of purely chiral oscillations, does 
not introduce relevant modifications to the flavor conversion formula. Otherwise, by taking 
the non-relativistic limit, with a momentum distribution sharply peaked around |pj <C 
the complete oscillation probability formula can be written as 

P(^, L - f) = ^ £ { ^ ( 1 " M 2 )aS '" ) - (cos [mPt] - cos [m^]) 2 } (48) 



s = 1,2 
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which, despite not being phenomenologically verifiable, introduces a completely different 
pattern for flavor/chiral oscillations. 

Meanwhile, a more interesting interpretation is provided when we analyze the second 
case, where the propagating momentum p is orthogonal to the magnetic field B. Now, 
in the flavor/chiral oscillation formula, the effects of the external magnetic field can be 
discriminated from the mass interference term, in the sense that we can write the complete 
expression for Dfco(£) as 

DFOO(t) = i| ^9 2 (P){(»8[|o">|(]co 8 [|al'>|(]) 

\ + " + flf ) c«i(4" - 4") t] + (i - p2+ f,f ) ™[(4" + 4") t] 



m (!) m ( 2 ) 

+ 



fc o fc o 



cos[( £ W +ei 2) ) t]cos[(|aW| - |a^|) t] 

- cos^ - 4 2) ) t] cos[(ja«| + |a( 2 )|) t]] } , (49) 

where we have used the correspondence of and with m''l Again, the ultra- relativistic 
limit reduces the impact of the modifications to residual effects which are difficultly de- 
tectable by experiments. By taking the non-relativistic limit and following the same proce- 
dure for obtaining the Eq. (48), the complete oscillation probability formula can be written 



as 



P(v a ,r - u p , L ;t) = (cos 2 [m^t] + cos 2 [m^t] 

-2cos[(|a (l) | - |a (2) |)t]cos[m (l) t]cos[m (2) t]) . (50) 

Turning back to the starting point, if we had postulated a wave-packet made up exclu- 
sively of positive frequency plane-wave solutions, the oscillation term which appear as a 
sum of mass-eigenstate energies would have vanished. The new oscillations have very high 
frequencies. Such a peculiar oscillating behavior is similar to the phenomenon referred to 
as zitterbewegung. In atomic physics, the electron exhibits this violent quantum fluctuation 
in the position and becomes sensitive to an effective potential which explains the Darwin 
term in the hydrogen atom [38]. It reinforces the argument that, in constructing Dirac wave 
packets, we cannot simply forget the contributions due to negative frequency components. 
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IV. CONCLUSIONS 



To quantify some subtle processes which accomplish the (standard) flavor oscillation 
phenomena [15] and emerge via chiral oscillations of propagating neutrinos non-minimally 
interacting with an external magnetic field B, we have reported about some recent results 
on the study of flavor oscillation with Dirac wave-packets [32]. By taking into account 
the spinorial form of neutrino wave functions and imposing an initial constraint where a 
pure negative- chiral flavor-eigenstate is created at t = 0, for a constant spinor w, we have 
calculated the contribution to the wave-packet propagation due to the interference between 
positive and negative frequency solutions of the Dirac equation with non-minimal coupling 
and, finally, we have obtained a complete expression for the oscillation probability. 

It effectively represents the formally accurate way for deriving the expression for the neu- 
trino spin-flipping in magnetic fields related to chiral oscillations and coupled with a flavor 
conversion mechanism. It is also to be noted that in the above discussion we have assumed 
that neutrinos are Dirac particles, thus making the positive- chiral component sterile. If the 
neutrino was a Majorana particle, it should not have a magnetic moment, obviating the 
spin-flipping via magnetic field interactions but still allowing the (vacuum) chiral conversion 
possibility via very rapid oscillations (zitterbewegung) . 

We have confirmed that the fermionic character of the particles modifies the standard 
oscillation probability which is previously obtained by implicitly assuming a scalar nature 
of the mass-eigenstates. At the same time, it allows us to correctly determine the origin and 
the influence of chiral oscillations and spin-flipping in the complete flavor conversion result. 
Strictly speaking, we have obtained the term of very high oscillation frequency depending on 
the sum of energies in the new oscillation probability formula which, in case of Dirac wave- 
packets, represents some modifications that introduce correction factors which are negligible 
in the UR limit. Our future perspectives concern with deriving the flavor coupled with chiral 
conversion expressions for neutrino moving in the background matter by supposing that the 
magnitude of some experimentally (implicitly) observable matter effects could be quantified 
(and eventually detected) in this Dirac wave-packet framework. In particular, under the 
point of view of a phenomenological analysis, we cannot discard the possibility of existing 
larger magnetic momentum (10~ 12 //s) for the electron/muon neutrinos deduced from some 
extensions of the minimal standard model [41, 42] which eventually can be implemented in 
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a subsequent study. 

We know, however, the necessity of a more sophisticated approach is understood. In fact, 
the derivation of the oscillation formula should resort to field-theoretical methods which, 
meanwhile, are not very popular. They are thought to be very complicated and the existing 
quantum field computations of the oscillation formula do not agree in all respects [16]. The 
Blasone and Vitiello (BV) model [9, 11] to neutrino/particle mixing and oscillations seems 
to be the most distinguished trying to this aim. But still with Dirac wave-packets, the flavor 
conversion formula can be reproduced [32] with the same mathematical structure as those 
obtained in the BV model [9, 11]. Moreover each new effect present in the oscillation formula 
can be separately quantified. 

Just to summarize, we would not have not been honest if we had ignored the complete 
analysis of the general case comprised by Eqs. (10-13) where we had not yet assumed an 
arbitrary spatial configuration for the magnetic field. Meanwhile, such a general case leads to 
the formal connection between quantum oscillation phenomena and a very different field. It 
concerns with the curious fact that the above complete (general) expressions for propagating 
wave-packets do not satisfy the standard dispersion relations like E 2 = m 2 + p 2 excepting by 
the two particular cases where E s (p) 2 = m 2 s +p 2 for p x B = or e'l = m 2 + p 2 for p- B = 0. 
In principle, it could represent an inconvenient obstacle which forbids the extension of these 
restrictive cases to a general one. However, we believe that it can also represent the starting 
point for discussing dispersion relations which can be incorporated into frameworks encoding 
the breakdown (or the validity) of Lorentz invariance. 
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APPENDIX A: FLAVOR OSCILLATIONS WITH WAVE-PACKETS - COMMON 
POINTS BETWEEN SCALAR AND FERMIONIC PRESCRIPTIONS 



By reviewing the flavor conversion process for which we denominate a scalar prescription 
[13], we notice that the time-evolution of flavor wave-packets can be described by 

4>(t,x) = (p 1 (t,x)cos6u 1 + <p 2 (t,x)sm9u 2 

= [</>!(*, x) cos 2 9 + 4> 2 (t,x) sin 2 0} u a + [^(t^x) — (f> 2 (t,x)] cos 6> sin # i/g 

= (f> a (t,x',e)u a + (f>p(t,x',e)i/ fi , (Ai) 

where v a and v$ are flavor-eigenstates and v x and u 2 are mass-eigenstates. The probability 
of finding a flavor state v$ at the instant t is equal to the integrated squared modulus of 
the Vp coefficient 



P(v a ^"f>;t) = J d 3 x \^(t,x;6)\ 2 = s -^{l-Fo(t)}, 



(A2) 



where Fo(t) represents the mass-eigenstate interference term given by 

Fo(t) = Re 



J d A X(j)\{t,x)(j) 2 {t,x) . (A3) 
As an illustrative example, we can consider gaussian wave-packets given at time t = by 

exp[ip 4 z], (A4) 



( 2 V 


x 2 ' 


— H ex P 




\na 2 J 


.~~a 2 '. 



where s = 1, 2. The wave functions which describe their time-evolution are 

dp 



(2tt) 3 



<f{p- Pi) ex P H #(P> mjt + ip- x], 



(A5) 



(p-pj 2 a 2 



&(t,a:) 

where E(p,m l ) = (p 2 + m 2 ) 3 and 

i 

y?(p — pj = (27ra 2 ) 4 exp 

In order to obtain the oscillation probability, we can calculate the interference term Fo(t) 
by solving the following integral 

dp 



■ <f{P - PiMp - P 2 ) ex P H a #(p) *] 



exp 



Ap) 5 



dp 

(2^ 



^ 2 (p-p )expHA£(p)t], (A6) 
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where we have changed the ^-integration into a p-integration and introduced the quantities 
A\p\ = p 1 — p 2 , P = \{Pi +P 2 ) an d AE(p) = E(p,m 1 ) — E(p,m 2 ). The oscillation term is 
bounded by the exponential function of a A\p\ at any instant of time. Under this condition 
we could never observe a pure flavor-eigenstate. Besides, oscillations are considerably sup- 
pressed if oA|p| > 1. A necessary condition to observe oscillations is that aA|p| <C 1. This 
constraint can also be expressed by S\p\ A|p| where 5\p\ is the momentum uncertainty of 
the particle. The overlap between the momentum distributions is indeed relevant only for 
5\p\ 3> A|p|. Consequently without loss of generality, we can assume 

Fo(t) = R e y^^(p- p ) exp [-i AE(p) t] } . (A7) 

In literature, this equation is often obtained by assuming two mass-eigenstate wave- 
packets described by the "same" momentum distribution centered around the average mo- 
mentum p = p () . This simplifying hypothesis also guarantees instantaneous creation of a 
pure flavor eigenstate v a at t — 0. In fact, for 0i(O, x) = 4> 2 (0, x) we get from Eq. (Al) 



(f> a {0,x,6) 



( 2 V 


X 2 ' 


exp 









exp [ip ■ x] (A8) 



and 0/3(0, x, 6) = 0. 

The time-evolution of a spin one-half particle follows an analogous prescription. Once we 
have introduced the fermionic character in the study of quantum oscillation phenomena, we 
can use the Dirac equation as the evolution equation for the mass-eigenstates represented 
by ifti(t, x) so that the natural extension of Eq. (A8) reads 

V> a (0, x, 6) = <p a {0,x, 6)w (A9) 

where we assume w is the constant spinor which satisfies the normalization condition w^w = 
1. 
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